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ABSTRACT 

Interacting galaxies favor the formation of star clusters but are also suspected to affect 
their evolution through an intense and rapidly varying tidal field. Treating this complex 
behaviour remains out-of-reach of (scmi-)analytical studies. By computing the tidal 
field from galactic models and including it into star-by-star iV-body simulations of star 
clusters, we monitor the structure and mass evolution of a population of clusters in 
a galaxy major merger, taking the Antennae galaxies (NGC 4038/39) as a prototype. 
On the long timescale (~ 10 9 yr), the merger only indirectly affects the evolution of 
clusters by modifying their orbits in or around the galaxies: the mass-loss of clusters 
in the merger remnant is faster, while clusters ejected in the tidal debris survive much 
longer, compared to in an isolated galaxy. The tidal perturbations of the galactic 
collisions themselves are too short lived and not strong enough to significantly influence 
the structure and dissolution of realistically dense/massive star clusters. 

Key words: open clusters and associations: general - galaxies: star clusters - galaxies: 
interactions - methods: numerical 



1 INTRODUCTION 

For 40 years n ow, since the pioneer work of 
iToomre fc Toomrd (|l972t ). many studies focussed on 
the role of mergers on galaxy e v olution (among many 



mergers on galaxy evolution (among many 
Elmegreen et al.l 1 19951 : iBournaud et~akl |200S| ; 



others, _ _ 

lLotz et all 120081 ). Hydrodyna mical simulations explore d 
the process of star formation (|Kravtsov fc Gnedinl [2005h 
with the goal of interpreting the observational data of 
(Ultra-) Luminous Infra-Red Galaxies, and formation of 
star clu sters often detec t ed sh o rtly after a major merge r 
event l| Ashman fc ZepJ 1 19921 : iBrodie fc Straderl |2006h . 
Numerical simulations confirmed indeed that galaxy-galaxy 
interactions could significantly increase the star cluster 
formation rate (e.g. |Pi Matteo et al.ll2007l ; iKarl et al.ll2010l : 
iTevssier et al.ll2010t ). However, at the stage of the conversion 
of gas into stellar material, a star-by-star exploration in 
a full galactic context, but also the subsequent dynamical 
evolution, remain technically out-of-reach. 

However, the fate of these clusters after their for- 
mation remains uncertain. This is because an evapo- 
rating/dissolving cluster in a galaxy beyond the local 
group leaves no clear observational signature for us to 
find. Numerical studies are challenged because they must 
account for a wide range of scales, from the internal 
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dynamics of clusters (i.e. binary stars and close encounters) 
to the scale of galactic interactions. This issue has been 
circumvented by simplifying either one or the other side of 
the problem. On the one hand, an idealized galaxy with 
an analytical description of the tides has been intensively 
used to improve our understanding of the dynamical 
evolu t ion of non-isolated clusters ( e.g. Vesperini fc H cggic 



1997; Baumgardt & Makino 2003; Penarrubia et al 



i Kiipper et all l20ld ; iBerentzen fc Athanassoula 



2009; 



20121 ; 



Iwebb et alj|2013l h" On the other, a semi-analytical descrip- 
tion of star clusters in a cosmological context allowed to 
probe the effect of complex tides on stellar p opulations (e.g. 
iPrieto fc Gnedinll2008l ; iKruiissen et al.ll201ll l. 



Recently, iRenaud. Gieles fc~ Boilvl (|201ll . hereafter 
RGB11) merged the two approaches by integrating any tidal 
effect into TV-body simulations of star clusters, allowing for 
the exploration of the evolution of star clusters within com- 
plex, time-dependent galactic potentials, like those of merg- 
ers. The goal of this Lettre is to apply this method to a 
large number of clusters in rapidly-varying tidal fields, and 
to derive the differences in their evolution and dissolution 
rate between isolated galaxies and mergers. Our approach 
focusses first on a few cases to better understand the phys- 
ical phenomenon of dissolution itself, and then on an entire 
population of stars clusters. This work does not intent to 
reproduce realistic clusters, but rather to shed light on the 
physical effect of complex time-evolving tides, leaving aside, 
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for the time being, other aspects like mass segregation and 
stellar evolution. 



2 METHODOLOGY 

2.1 Numerical method 

The evolution of our cluster models is com puted wi t h the 
codjB MB0DY6tt, w i dely based on NB0DY6 jAarsethl |2003l ; 
iNitadori fc Aarsethl l2012h and presented in iRGBlli and 
briefly summarised here. First, a simulation of a galaxy is 
performed: each particle represents a possible star cluster. 
(This step can be done using any type of gravitational code, 
e.g. tree-code.) One particle is followed alon g its orbit and 
the a ssociated tidal tensors are extracted (see Renaud et al.l 
2008), providing a complete space and time description of 
the tides for this cluster candidate. Next, in a second sim- 
ulation, NB0DY6tt creates a star-by-star iV-body model of 
a star cluster, reads the tensors, interpolates them in time, 
computes the tidal forces at the positions of the -/V stars 
of the cluster and adds them to the internal gravitational 
force due to the N — 1 other stars. A significant fraction of 
the computation is speed-up thanks to the use of Graphics 
Process ing Units (GPUs) . Details and tests can be found in 
IRGBlli . In this Lettre, we repeat this exercise by considering 
various orbits of clusters, i.e. various tidal fields, in a merger 
and in an isolated galaxy. To focus on the tidal effect and on 
the dissolution of the cluster population, we do not include 
any prescription on the formation of our clusters. In total, 
724 iV-body simulations have been run, on the GPU cluster 
Curie hosted at the Tres Grand Centre de Calcul (TGCC). 

2.2 Galaxies, orbits and clusters 

We computed the tidal tensors using a tree-code 
along several orbits in a purely gravitational model of 
N GC 4038. (The met hod and simulation details are given 
in lRenaud et aLlbOOSl .) This galaxy has been modeled both 
in isolation and with its equal-mass companion (NGC 4039) , 
the pair reproducing the morphology and kinematics of the 
Antennae galaxies. This allows a direct comparison between 
the isolated galaxy and the major merger. 

As the galactic collisions proceed, the tidal history gets 
punctuated by "events" corresponding to the impact of the 
two galaxies. These events however are independent of the 
intrinsic ev olution of the clusters, itself pac ed by two-body 
relaxation (|Henonl Il96ll ; iGieles et al] |2010l. To avoid any 
artificial synchronisation between the galaxy and the clus- 
ter, we start the integration of the cluster models at different 
times with respect to the collisions, i.e. at different positions 
along the orbits. This way, our sample comprises an entire 
population of young, intermediate and old clusters. How- 
ever, all the simulations are stopped 5 Gyr after the first 
galactic pericenter passage: for longer timelapses, cosmolog- 
ical aspects stepping in the galaxy evolution should be taken 
into account, which is out of the scope of the present study. 
With 32 orbits and five "birth" epochs along them, we can 
monitor histories representative of all cases observed in a 
merger: clusters confined to the galaxy nucleus, or ejected 

1 http:/ /irfu. cea.fr/Pisp/florent.renaud/nbody6tt.php 



in the close vicinity of the central region, in tidal bridges or 
tails, falling back and orbiting the remnant, or being defini- 
tively ejected in the intergalactic medium, and this at differ- 
ent stages of their evolution. Half of the 32 orbits have been 
selected to be confined to the galactic disc (before the inter- 
action occurs) while the others reach positions in the stellar 
halo above and below the plane of the disc, on very inclined 
orbits. The maximum apocenter corresponds to the radius 
of the stellar disc (20 kpc) such that our survey represents 
a full range of orbital energies (see Section \3. 2 \ . 

For simplicity, in each cluster, all stars have been given 
the same mass (1 Mq), to avoid mass segregation and stel- 
lar evolution effects and to focus on the tides. The initial 
positi ons and velocitie s of the stars were based on a viri- 
alised iPlummerl (|l91ll ) model. On each orbit and for each 
birth epoch, we set up clusters with different masses and 
densities obtained either by changing the initial number No 
of stars per cluster from 4000 to 32000 for a virial radius 
r v = 1 pc, or by changing r v (down to 0.5 pc) for a given 
Nq. The half-mass densities vary from 10 3 to 10 4 Mq pc~ 3 , 
which makes our models comparable in mass and density to 
Westerlu nd 1, NGC 3603 or even th e Arches for the dens- 
est ones (jPortegies Zwart et al.l [20101 ) . However, because of 
early evolution, the densities get lower (down to two orders 
of magnitude in the extreme case) at the time of the galactic 
interaction. 



3 RESULTS 

By looking at the mass evolution curves of all individual 
clusters in our survey, we found shows that every case is 
unique. However, it is possible to extract families of clusters 
sharing a similar tidal history. 



3.1 Individual cases 

First, in Fig. [I] we show two orbits from our survey, in order 
to analyse typical situations that may occur at galactic scale 
and their implications for the clusters. Only the simulations 
with a birth epoch of ~ 1.5 Gyr before the first pericentre 
passage are shown here, for the sake of clarity, but the other 
runs give comparable results. 

On the top row case (panels abc.l), the orbital eccen- 
tricity is rapidly increased from ~ 0.05 to ~ 0.30 by the 
galactic collision. However the trajectory remains closed and 
fairly regular. The resulting mass evolution shows the typ- 
ical "staircase" decline, the mass-loss becoming more rapid 
at small galactocentric radii because of strong gravitational 
tides and a significant centrifugal force. The amplitude of 
this periodically accelerated mass-loss increases with low- 
mass clusters. This is also found fo r eccentric orbits (e.g. 
iBaumgardt fc Makinrj I2003L [rGBIII ). Despite these short- 
period (~ 10 8 yr) changes induced by the merger, the sec- 
ular evolution of this cluster does not strongly differ from 
that it has in an isolated disc (dashed lines on Fig. [T]). In 
that sense, the merger does not affect the mass-loss rate 
of clusters on this particular orbit. In Section T3.2I we show 
that this is not a representative example, in fact, because the 
mass-loss rate of clusters that remain bound to the merger 
remnant increases on average. 

On the contrary, the orbit of our second example (panels 
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t [Gyr] x [kpc] 

Figure 1. Evolution of clusters along two orbits arbitrarily se- 
lected. Panels (a): normalized mass of the clusters, without tides 
(dotted curves), in the isolated galaxy (dashed) and in the merger 
(solid). The colour indicates the initial mass of the cluster (4000, 
8000, 16000 and 32000 stars for blue, green, red and black respec- 
tively, all with an initial virial radius of 1 pc). Vertical dotted lines 
mark the galactic pericentre passages; the second one corresponds 
to the final coalescence. Panels (b): Galactocentric distance. Pan- 
els (c): Orbit of the clusters in the orbital plane of the merger. 
(The clusters start from the top-right corner. A diamond and a 
triangle indicate the positi on of the clusters a t the times of the 
two galactic collisions. See iRenaud et"ai1l2009l for details on the 
galactic simulation.) 



abc.2) is much more altered by the galactic collision. At the 
first galactic pericentre passage (t = 0) , the orbit follows the 
creation and expansion of a galactic tidal tail and does not 
fall back into the nucleus. Being in a shallow potential, the 
clusters on this trajectory experience a very weak tidal field, 
and thus a slow mass-loss (the mass-loss rate for t > 0.3 Gyr 
is comparable to that of an isolated cluster, i.e. most escape 
happens because of close encounters in the core rather than 
a relaxation driven evaporation over the tidal boundary), 
whereas in the isolated disc, the orbit remains bound and the 
clusters dissolve in a few xlO 9 yr. In this case, the merger 
protects the clusters by sending them to a zone of weaker 
tides. 



3.2 Cluster population 

To generalise these results, we now consider all our clusters, 
but tell apart two sub-populations: the clusters always stay- 
ing within 20 kpc from the center of the galaxy (i.e. about 
the radius of the isolated disc) and the others, ejected to 
the tidal debris. About 30% of the orbits fall in this sec- 
ond category, which corresponds to the stellar mass fraction 
lost by the progenitor galaxies of the Antennae during the 
interaction. This demonstrates that our sample of orbits is 
statistically representative. (For the isolated disc, this dis- 



Figure 2. Evolution of the normalized mass averaged by telling 
apart the clusters ejected to large galactocentric radii by the 
galaxy-galaxy collision (solid thick line) from those remaining 
close to the center (solid thin line). The dashed lines represent 
all the clusters in the isolated disc and the dotted lines denotes 
the mass-loss without tides. As in Fig.[T] the colour indicates the 
initial mass of the cluster. 



tinction is not made since all the clusters remain bound to 
their host galaxy.) 



3.2.1 Mass-loss 

Fig. [2] shows the average mass for the two groups and for 
two birth epochs. The "staircase" behaviour is smoothed 
out because all the orbits in our sample which might ex- 
hibit it are not synchronised. As expected from the exam- 
ple of the previous Section, the ejected clusters (thick solid 
line) live much longer as the tidal field they experience is 
weaker than in the merger remnant and the isolated galaxy. 
The knees in their mass evolution, marking a change of tidal 
regime, occur at first pericentre passage, when the tidal tails 
are created. However, at that time, the non-ejected clus- 
ters only experience the effect of the cumulated potentials 
of the two galaxies (before they separate and fall-back onto 
each other), which is too short-lived to significantly affect 
the mass- loss (see Section T3.3p . However, at the final coales- 
cence (t = 0.3 Gyr) the curve of the mass evolution yields 
a sharp knee. Contrarily to the particular case presented in 
the top row of Fig. [1] on average the mass of the non-ejected 
clusters decreases much more rapidly in the merger than in 
the isolated disc. Clusters in the central region of a merger 
remnant are very likely to dissolve faster (within a few Gyr, 
for the masses considered here) than in an isolated galaxy, 
partly because of an increased average orbital eccentricity 
from 0.09 ±0.03 before the first collision, to 0.39 ±0.12 after 
the merger. 

Finally, Fig. [3] compares the evolution of the mass of 
the entire population of a galaxy, independently of the clus- 
ter density, mass or birth epoch, but still telling apart the 
ejected clusters (tidal debris) from those staying in the rem- 
nant. The pre-merger history (t < 0) is not shown here 
because it encompasses the artificial creation of clusters 
at arbitrary times and thus, is not physically relevant for 
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Figure 3. Dissolution of the entire population of clusters (all 
densities, all birth epochs and all orbits) for our four families 
of tidal histories since the first galactic collision. The mass is 
normalized to its value at the first pericenter passage. The number 
in parenthesis is the characteristic timescale r of the exponential 
decay of the mass (M ex e - */ 1 "). 



Figure 4. Evolution of the half-mass radius of our clusters. The 
curves are stopped when the corresponding average normalized 
mass is less than 0.15, making the half-mass radius very Poisson- 
noisy. Colors and linestyles are as in Fig. [2] 



our comparisons. In all cases, the average mass decays as 
M oc e _t//T up to a very late stage of the evolution. The 
dissolution rate — 1/r of the population in the tidal debris is 
comparable (yet a bit higher) to that of the tide- free clusters 
and ~ 4 times lower than in the merger remnant. The iso- 
lated galaxy exhibits an intermediate dissolution rate, twice 
lower than that of the merger remnant. From the results 
in Fig. El we thus conclude that the average mass-loss rate 
increases for clusters that remain bound to the merger rem- 
nant. For some clusters the mass-loss stays comparable (e.g. 
the example in the top panel of Fig. [TJ, but it increases for 
a majority of clusters. 



3.2.2 Structural evolution 

In this Section, we focus on the structural cha nges along 
the evolution of our models. iGieles et al.l (|201ll ) explained 
that, after a relaxation-driven expansion phase, clusters in a 
constant tidal field experience a contraction phase. In Fig. [4] 
we show that clusters in time- varying fields also experience a 
comparable two-phase evolution. The transition from expan- 
sion to contraction occurs half-way th e lifetime of the clus- 
ter, as predicted bv lGieles et al.l l|201ll) for steady tides. This 
demonstrates that the variations of the tidal field (~ 10 s yr), 
i.e. rapid compared to the lifetime of the cluster (~ 10 9 yr), 
do not modify qualitatively the behaviour from what is ex- 
pected in constant tides. 

As noted above, the clusters remaining in the central 
region of the merger (thin solid lines) experience stronger 
tides than in the isolated galaxy (dashed lines) which makes 
their contraction more important. The clusters ejected into 
the tidal debris (thick solid lines) follow a quasi tide-free 
evolution and thus miss the tidally-driven contraction phase. 
In Fig. [5] we plot the evolution of the half-mass radius of 
one cluster set on the orbit shown in Fig. \T\ (bottom panels). 
At the time of the first collision, the expansion phase of this 
cluster is over and the contraction has just begun. However, 
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Figure 5. Evolution of the half-mass radius of a cluster sent in 
the tidal debris (one cluster of the group shown by the blue solid 
thick line in Fig. [4}. The new tidal regime starting at t = allows 
for a resumption of the expansion phase and leads to an extended 
cluster. 



the sudden change of tidal regime due to the ejection to the 
tidal tails stops the contraction and allows for a much longer 
expansion. 

3.3 Role of collisions 

In the previous Section, we showed the indirect changes im- 
plied by the merger on its population through orbital mix- 
ing, i.e. a long-term effect. However, the tides are expected 
to have their maximum influence at the time of the galac- 
tic pericenter passages, when the gravitational potentials of 
the still undisturbed galaxies most overlap. At the collision 
(t — 0), the velocity difference between the two galaxies 
is ~ 200 km s" 1 , therefore a significant fraction of their 
masses overlaps for a comparable or shorter timelapse than 
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the crossing times of our clusters (~ 10 — 100 Myr). Further- 
more, it should be noted that the geometry of the mergers 
does not necessarily account for a global strengthening of 
the tides: one galaxy can locally compensates for the tides 
of the other by flattening the gravitational potential]. In 
short, during the collision, the tidal field could be locally 
weaker than expected, but above all, the combination of the 
effects of the two galaxies is short-lived compared to the 
crossing time of their clusters. 

In our survey, the more fragile clusters loose only ~ 
2% of their current mass during the first galactic colli- 
sion. For instance, the clusters like that of Fig. [S] have an 
initial half- mass density of 1000 Mq pc -3 , which becomes 
~ 10 Mq pc~ 3 at the time of the collision. Even on such 
low density objects, the direct effect of the collision is very 
mild. It is visible in the mass-loss curves (see the blue curves 
in Fig. [T]a.l and Fig. [3] around t = 0) but is negligible with 
respect to the classical rapid mass-loss observed when the 
cluster reaches the pericenter of its orbit (i.e. the steep part 
of the "staircase" shape mentioned above). Although the 
orbital changes occurring in the merger remnant makes a 
similar analysis much more difficult at the second collision 
(t ~ 0.3 Gyr), we expect qualitatively this behaviour to re- 
peat itself. In fine, the galactic collisions themselves have no 
direct impact on the mass-loss of clusters. 



4 DISCUSSION AND CONCLUSION 

We have studied the tidal effect of a galaxy major merger 
on its population of star clusters. Our main findings are: 

• the mass-loss rate of the clusters that remain bound to 
the merger remnant is higher than before the final coales- 
cence phase. The clusters populating the tidal debris survive 
much longer, similarly to tide-free cases. 

• the two-phases evolution of expansion and contraction 
of a cluster also exists in complex, time-varying tidal fields. 
The structural evolution does not strongly deviate from that 
of a constant tides case. 

• at the time of the galactic collisions, the tides are too 
weak (although they reach their maximum intensity) and 
too short-lived (< 10 7-8 yr) to have a significant influ- 
ence on the clusters. However, by modifying the orbit of 
the clusters, they indirectly affect their mass-loss over long 
timescales. 

By focussing on tidal effects, we have not considered 
stellar evolution, which influences the evolution of the mass 
and the structure of clusters. Stellar evolution causes clus- 
ters to expand faster in the first few 100 Myr, making clus- 
ters less dense and possibly more vulnerable to tidal forces. 
However, because in this study we have considered clusters 
with different physical ages, some of which already reached 
the tidal density before the merger, we do not expect that 
our overall conclusions will be affected by the inclusion of 
stellar evolution. 

2 In the idealised case of two point-mass galaxies, the effective 
potential has a saddle point where the tidal forces balance the 
gravitation of the cluster. For more complex galaxies, this situa- 
tion may occur over larger volumes, depending on the shape of 
the potential and the geometry of the encounter. 



Furthermore, we have not taken the formation of the 
clusters into account. High resoluti on hydrodynamical sim- 
ulations of the same galaxy model jTevssier et al.ll2010T l re- 
veal that the star formation activity increases by a factor up 
to ~ 20 shortly after the pericentre passages when the tur- 
bulent interstellar medium is properly resolved. However, 
galaxy-scale simulations still treat star formation through 
subgrid recipes and most of them lack the resolution to 
properly describe the initial state of clusters. For theses rea- 
sons, our approach decouples the dissolution of the clusters 
from their formation, considering separately the artifacts 
from both sides. Hence, our results do not represent the 
age or mass functions of clusters, but their mass-loss rate. 
Studying the formation sites of the objects presented here 
(~ 10 3-4 Mq clusters) in a fully consistent galactic context 
by resolving star forming cores and the stellar feedback down 
to ~ 1 pc is our next step (Renaud et al., in prep). 
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